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Thermal statistical models are simple and effective tool to describe particle production in high 
energy heavy ion collision. It is shown that for higher moments finite volume corrections become 
£ — | important observable quantities. They make possible to differentiate between different statistical 

f^*) , ensembles even in the thermodynamic limit. 

o 

(N 



o 



I. INTRODUCTION 



Particle production yields are nicely reproduced by thermal models, based on the assumption of 
(N ■ nonintcracting gas of hadronic resonances [l|. This simplicity can be however misleading as particle 

yields and particle ratios are not very sensitive to the underlying model. The main ingredient of 
the statistical models are probability densities, which allow to extract the whole physical informa- 
• tion. The only way to reproduce those probability distributions is by means of higher and higher 

' probability moments. These moments are in fact the only quantities which are phenomenologically 

available and can be used for the verification of theoretical predictions. 

Particle yields in heavy ion collision are the first moments, so they lead to rather crude com- 
parisons with the model. Fluctuations and correlations are second moments so they allow for the 
better understanding of physical processes in the thermal equilibrium. 

Fluctuations and correlations measured in heavy ion collision processes give better insight into 
dynamical and kinematical properties of the dense hadronic medium created in ultrarelativistic 
heavy ion collisions. Systems under considerations are in fact so close to the thermodynamic limit 
that final volume effects seem to be unimportant — at least when productions yields are considered. 

The aim of the paper is to show that finite volume effects become more and more important 
when higher moments, e.g. correlations and fluctuations are considered. 

A preliminary analysis of the increasing volume effects was given in @, Q ■ It has been rigorously 
shown an influence of 0(1/V) terms for a new class physical observables — semi-intensive quan- 
tities @. Those results completely explained also ambiguities noted in related to "spurious 
non-equivalence" of different statistical ensembles used in the description of heavy ion collision 
processes. 

We start from the simple example of the standard statistical physics. Is is shown that even in 
that case a notion of semi-intensive quantities is relevant for the physical situation. In the next 
step we consider an abelian symmetry corresponding to one conserved charge. 



II. CHOICE OF VARIABLES 



In the thermodynamical limit the relevant probabilities distributions are those related to densi- 
ties. These distributions are expressed by moments calculated for densities — not for particles. In 
the practice, however, we measure particles — not densities as we do not know related volumes. 
Fortunately, volumes can be omitted by taking corresponding ratios. 

Let us consider e.g. the density variance An 2 . This can be written as 
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By taking the relative variance 

An 2 _ (N 2 ) - (N) 2 
(n)* ~ {N)~ 2 ' 

volume-dependence vanishes. 

A. Semi-intensive variables 

A special care should be taken for calculations of ratios of particles moments. Although moments 
are extensive variables their ratios can be finite in the thermodynamic limit. These ratios are 
examples of semi-intensive variables. They are finite in the thermodynamic limit but those limits 
depend on volume terms in density probability distributions. One can say that semi-intensive 
variables "keep memory" where the thermodynamic limit is realized from. 

Let consider as an example the scaled particle variance 

{N 2 )-{N) 2 _ (n 2 )-{n) 2 
(N) (n) ■ 

The term 

in 2 ) (n) 2 
(n) • 

tends to zero in the thermodynamic limit as 0(V^ r ). So a behavior of the scaled particle variance 
depends on the 0(V _1 ) term in the scaled density variance. A more detailed analysis of semi- 
intensive variables is given in Q. 

To clarify this approach let us consider a well known classical problem of Poisson distribution 
but taken in the thermodynamic limit. 



III. GRAND CANONICAL AND CANONICAL ENSEMBLES 

A. Poisson distribution in the thermodynamic limit 

Let us consider the grand canonical ensemble of noninteracting gas. A corresponding statistical 
operator is 

D=— -. r- (1) 

This leads to the partition function 

Z(V,T, 7 )= e ze \ (2) 
where z is one-particle partition function 

Z(T ' V) = J^J ^ = VZa{T) ' (3) 

A 7 parameter (= /3^) is such to provide the given value of the average particle number 
(N) = V(n) . This means that 



Particle moments can be written as 



e 7 = XZL . ( 4 ) 



1 r) k 7 
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The parameter 7 is taken in final formulae as a function 7((n), Zq) from Eq ([¥]). 
The resulting probability distribution to obtain N particles under condition that the average 
number of particles is (N) is equal to Poisson distribution 



We introduce corresponding probability distribution V for the particle number density n — N/V 



V {n) (n;V) = VP v{n) (Vn) =V 



(V(n)) Vn 



-V(n) 



T(Vn+l) 

For large Vn we are using an asymptotic form of Gamma function 



(G) 



T(Vn+l) ~ V2ir(V 



iVn-l/2 e -Vn 



1 



1 



l2Vn 



0(V~ 



This gives 



P<n,(n;V)~W» 1 (M) 
V2irn \ n J 



Vn 



1 



\2Vn 



0{V~ 



(7) 



This expression in singular in the V — > 00 limit. To estimate a large volume behavior of the 
probability distribution © one should take into account a generalized function limit. So we are 
going to calculate an expression 



(G) v = j dnG(n)V (n) (n;V) , 

where "P{ n }(n; V) is replaced by the asymptotic form from Eq ((7|). In the next to leading order 
in 1/V one should calculate 



V 1 ' 2 - 



1 



dn ^ e VS(n) _y-l/2 1 



(8) 



/2tt J n 1 / 2 12V2W n 3 / 2 

where 

S(n) — n ln(n) — n In n + n — (n) . 
An asymptotic expansion of the function (G) y is given by the classical Watson-Laplace theorem 
Theorem 1 Let I — [a, b] be the finite interval such that 

1. max5(i) is reached in the single point x = x , a < xo < b. 

2. f(x),S(x)eC(I). 

3. f(x),S(x) £ C°° in the vicinity ofxo, and S (xo) 7^ 0. 
Then, for A — > 00, A £ 5 C , there is an asymptotic expansion 



F[X] ~ e^^^CfeA-^ 1 



/2 



r(fc + 1/2) ( d_ 



2 k 



f(x) 



S(x ) - S(x) 
(x - x ) 2 



-fc-1/2' 



(9a) 
(9b) 



S e is here a segment | argz| ^ \ — e < -| in the complex z-plane. 
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To obtain 0(1/V) formula the first term in ([5]) should be calculated till the next to leading order 
term in X/V. For the second term it is enough to perform calculations in the leading order only. 
The first term gives the contribution 

v " 2 i / *^' r ™ = G <«>+ i^(w>+f£ G "(«>. (10f " 



and the second term gives 



So we have eventually 



V-W -±= / dn^§ e™<»> = j^G{{n)) , (10b) 
12V27T J n 6 ' 2 - 12(n)V 



(G) v = G((n)) + ^G ((«)) + 0(F- 2 ) , (11) 



for any function G. 

This gives us the exact expression for the density distribution ^ in the large volume limit 

V {n) {n- V) ~ S(n - (»)) + M *"(„ - („)) + 0(F" 2 ) . (12) 
We are now able to obtain arbitrary density moments up to 0(V~ 2 ) terms. 

(n k )v = J dnn k V {n) (n; V) = (n) k + fcil ( n )k-i + {V~ 2 ) . (13) 
We have for the second moment (intensive variable!) 



<- 2 }v = <-> 2 + ^ + <^- 2 ). 

This means 

An 2 = M^ . (14) 

as expected in the thermodynamic limit. 

The particle number and its density are fixed in the canonical ensemble so corresponding vari- 
ances are always equal to zero. The result (|14p can be seen as an example of the equivalence of 
the canonical and grand canonical distribution in the thermodynamic limit. This equivalence is 
clearly visible from the Eq (JT5J) where the delta function in the first term can be considered as the 
particle number density distribution in the canonical ensemble. 

A more involved situation appears for particle number moments (extensive variable!). Eq (|13p 
translated to the particle number gives 

(N k ) = V k (n) k + V k - X k< " k 2 ^ (n) fc ~ 1 + 0(V k - 2 ) , (15) 

One gets for the scaled variance (semi-intensive variable!) 

AN 2 . v 

W = 1 ' (16) 

what should be compared with zero obtained for the canonical distribution. 

The mechanism for such a seemingly unexpected behavior is quite obvious. The grand canonical 
and the canonical density probability distributions tend to the same thermodynamic limit. There 
arc different however for any finite volume. Semi-intensive variables depend on coefficients at those 
finite volume terms so they are different also in the thermodynamic limit. 
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B. Energy distribution 



It is interesting to perform similar calculation for the energy distribution in both ensembles. 
Energy moments and an average energy density can be written as 



(E k ) = (-1)' 



1 d k Z 
Z~dfW 



dz 
'dp 



One gets from Eq (fTTf 



k,Ak , T.fc-i^O- 1 )/ \fc-2(") d 2 z 



0(V k 



k-2\ 



The grand canonical energy density distribution follows 

P(e\{n), (e» = 6 (e - (e)) + M f ifT) *"( e - (e)) + 0{V- 2 ) 



2V 



K bo is given here as 



(n 



1 d 2 z 



z dp 2 



0=0((e)/(n)) 



For the canonical distribution a corresponding statistical operator is 



D = 



This leads to the partition function 



z(y,T) 



Nl Nl 

Internal energy moments are given by Eq (|17j) . In particular 



n dzo 
z Q dp 



For the energy moments one gets now 



(E k ) = V k (e) k + y fe -!^ l l( t ) k - 2 ~ 9 f ldz ° 



dp{- W ]+o{v 



rk-2\ 



A corresponding probability distribution is 



P(e|n, (6» = 5 (e - (e)) + — T& ^ )S (e - <e)) + 0(V 



where 72. is given here as 



2V 



n I dp \zq dp 



(17) 



(18) 



(19) 



(20) 



(21) 



(22) 



(23) 



(24) 



IV. HIGH ENERGY STATISTICAL PHYSICS 



Although the spirit and the philosophy of the statistical approach remains the same, ingredients 
of statistical models used in high energy problems are different. The main difference is that 
a number of particles is not longer conserved so we have no chemical potentials related to that 
quantity. The only nontrivial chemical potentials are those related to conserved charges, so the role 
of internal symmetries is a crucial one. For abelian charges, as electric charge or baryonic charge 
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an introduction of corresponding potential is rather obvious - it is just a Lagrange multiplier at 
a generator of the U(l) symmetry. For non-abelian symmetries, as e.g. for the isotopic SU(2) 
symmetry, problem is more involved. It appears Q that in such a case the only relevant chemical 
potentials are those related to so called Cartan subgroup - a maximal abelian subgroup of the 
given non-abelian group of internal symmetry. 

So for the hot hadronic gas a well approximated internal symmetry is SU (3) or SU (4) with a 
charm taken into account. On the top of it we have the exact 17(1) baryon number conservation. 
SU{3) symmetry leads to two chemical potentials, related to the third isospin component and to 
the hypercharge. Those two chemical potentials are supplemented by the baryonic potential. One 
uses for practical reasons another set of chemical potentials, which are linear combinations of the 
basic set. These are the electric charge and the strangeness chemical potentials together with the 
unchanged baryonic potential. 

Then, for the simplest case of an ideal hadron gas in thermal and chemical equilibrium, which 
consists of I species of particles, energy density e, baryon number density fig, strangeness density 
ns and electric charge density n/Q) read (H = c = 1 always) 



--gS"- +i >/* M $: }+B - ™ 

i 00 

" s =^g>' +i) / V{$} + , - 

i 00 

i=i o cxp I rp > + gi 

where Ei = (mf + p 2 ) 1 / 2 and mi, Bi, Si, fii, s, and gi are the mass, baryon number, strangeness, 
chemical potential, spin and a statistical factor of specie i respectively (we treat an antiparticle 
as a different specie). And ^ = BifiB + SiHs + Qi n Q, where /is, /is, and [Xq are overall baryon 
number and strangeness chemical potentials respectively. 

To get particle yields one should consider also entropy density s 



1 ' 7 p 4(E i -» i )ex V {*^} 



i=l 



To obtain the time dependence of temperature and baryon number and strangeness chemical 
potentials one has to solve numerically equations (|25|) - (|26p with s, ng, ng, and rig given as time 
dependent quantities. For s(t), n_e(£), and nQ(t) one obtains expressions form hydrodynamical 
calculations and ns = since we put the overall strangeness equal to zero during all the evolution. 

These equations, enriched by unstable particles effects, form a basic for successful calculations [l[ 
of relativistic heavy ion production processes concerning particle yields and rates. All calculated 
observables are here the first moments of related probability distributions. If we are going to get 
correlations and fluctuations predictions, we have to calculate second moments. This gives quite 
new effects, statistical ensemble dependent, as was shown in previous sections devoted to standard 
statistical physics approach. 



A. Statistical ensembles of high energy physics 



To make our considerations the simplest possible we consider the statistical model of a non- 
interacting gas constrained by the conservation of the abelian charge Q. The thermodynamic 
system of volume V and temperature T is considered to be composed of charged particles and 
their antiparticles carrying charge ±1 respectively. The requirement of charge conservation in the 
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system is imposed on the grand canonical or canonical level. The canonical level means the global 
charge conservation while the grand canonical level means the charge conservation on the average. 
The partition function of the above canonical and grand canonical statistical system is found to be 



2g(y,T) = Tr Q e-^= £ i_ rr =/ g (2«), (27a) 

N + — N_=Q ' +' 

Z GC {V,T) = Tre-^-^ = exp (2zcosh^) . (27b) 
where z is the sum over all one-particle partition functions 

*(T)= / = l~ 2 T^ 9i m*K 2 (^ ) =V«P), ^ 

and gi is the spin degeneracy factor. The sum is taken over all charged particles and resonances 
of mass rrii carrying the charge ±1. The functions Iq and K% are modified Bessel functions. The 
chemical potential fj, determines the average charge in the grand canonical ensemble 

This allows to eliminate the chemical potential from further formulae for the grand canonical 
probabilities distributions 



| = arcsin hf = m^ + ^f W . (29) 

In the canonical ensemble we have a system of volume V and total charge Q. In the grand 
canonical ensemble we have a system with volume V and average charge (Q) . Number of particles 
is not conserved in both ensembles. Number N_ of negative charged particles shall be extracted 
from the relevant probability distributions. The probability distribution Vq(N^,V) to have iV_ 
negatively and N + = N_ + Q positively charged particles is obtained @, [f| from the partition 
function (|27ap as 

2N-+Q i 

^ = + q)l /„(*) ■ (30) 

On the other hand in the GC ensemble with volume V and average charge (Q) the probability dis- 
tribution VPqAN-, Q,V) to find a system with a given charge Q and a given number of negatively 
charged particles N_ is expressed [2| as the product 

Pffi (N-,Q, V) = T$(N-,V) Vffi (Q, V) , (31) 

of the canonical particle number distribution Vq(N^, V) from Eq. (|30|) and the grand canonical 
probability distribution 



P<Q)(Q>tO=/«(2z) 



(Q) + V(Q) 2 + ^ 2 



2z 

to find the total charge Q in the system with the average charge (Q) 



-V«> 2 +4* 2 (32) 



B. The thermodynamic limit 

The thermodynamic limit is understood as a limit V — > oo such that densities of the system 
remain constant. So we have for the canonical ensemble 
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and 

in\ at ^ i \ N - 

(Q) oo, iV_ -> oo ; — = {q); — = n_ 

for the grand canonical ensemble. 

To formulate correctly the thermodynamic limit of quantities involving densities, one defines the 
following probabilities 

V C q {n-,V) := VP$ q (Vn-,V) , (33a) 

Pffi{n-,q,V) := V 2 V§fc (Vn_ , Vq, V ) , (33b) 

Pf 9 f(g,V) := W£&(Vg„V). (33c) 

We are going to proceed now in a similar way as in the former section. In a large volume limit 
one gets 

? C q {n^V) = P~(n_) + ii£(n_) + 0(V- 2 ) , (34) 
Pf q f(n-,q,V) = ^ ) (n_,g) + ii?f 9 f(n_, (? ) + 0(y- 2 ), (35a) 



Pf,) (9, = (<?) + (g) + 0(V~') . (35b) 

All functional coefficients can be obtained here using Laplace- Watson theorem. From the careful 
analysis one gets Q for the probability distribution (l34l) of the canonical ensemble 



P°(n-;V) = 8 (n_ - (n_)oo) + ^ q 2 + ± z 2 S ' (»- - (»-)<») 

+ I_|L_ «J" ( n _ - (n_)oo) + 0(1/V 2 ) : 

and for the probabilities distributions (|35[) of the grand canonical distribution 
P? q ?(q,n.;V) = 5 (n_ - (n_)«,) % - (?)) 



(36) 



7 00 «J" (n_ - (n_U 5(g - <g)) + 0(1/T/ 2 ) , 



2V 



(37a) 



Pf g f (g, VO = « (« - (g» + VW 2T 7 <s (g-(g)) + G(i/^ 2 ) , (37b) 



Pgf (n_, V) = * (n_ - (»»_)«,) + ^P<*" (n_ - (n_)oo) + G(l/F 2 ) . (37c) 
An average limiting density of charged particles 



^q 2 + Azl±q 

(n±)oo = = (38) 



is used in above formulae. 



C. Particle moments 



Probability distributions l[3"6"j) and (pJTl) allow to write compact expressions for for particle and 
charge distribution density moments of any order up to C(l/T^ 2 ) terms. For particle moments one 
gets 
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for the canonical ensemble moments and 

{ nif^{n ± )l + ^^(n ± )^. (40) 

for the grand canonical ensemble moments. 

Although those moments are density moments they can be expressed by directly observable 
variables. Using Eq ([38]) one gets moments as functions of q/zo ratio. This ratio is observable as 
it can be written as a function of the ratio of charged particles 

q 2 _ (^+)oq (N-)oo 

zl (N^oo (7V+) 00 



D. Semi-intensive variables 



Now we are in position to create some semi-intensive variables. They are finite in T-limit and 
have different values dependently on how the charge conservation is implemented in the description 
of the system. There is actually a broad class of variables. We take an an example 

_ (N k ) (N) k 

Sk - <jv)*-i (41) 

Indeed from (f39|) and ([37]) one gets canonical and grand canonical values for positive(negative) 
particles in the thermodynamic limit (denoted as T-limit in the subsequent formulae) 



T- Bm5f = k{k - l) ^ 2 ±^1 , (42a) 



for the canonical ensemble, while in the grand canonical ensemble 

T-limS^fcl). (42b ) 

The scaled variance is just a special case of Sk corresponding to k = 2. 

Another examples are classes of variables closely related to cumulant or factorial cumulant 
moments [3J or susceptibility ratios. Let define p— th order susceptibility 



d p \nZ 
dfjP 



One can easily check that the ratios 



IC p -,r = ^ , (43) 



are semi-intensive quantities. 

One can also construct more involved semi-inclusive variables having a finite T-limit behavior 
which are determined by higher order asymptotic terms of the corresponding probability distribu- 
tions. 



V. CONCLUSION 



We have discussed the differences in the asymptotic properties of the probability functions for a 
system with an exact, that is canonical, and with an average, that is grand canonical, implemen- 
tation of charge conservation. We have shown that in the thermodynamic limit the corresponding 
probability distributions in the grand canonical and canonical ensembles coincide and are described 
as generalized functions. This property is a direct consequence of the grand canonical and canoni- 
cal ensemble equivalence in the thermodynamic limit. However, the first finite volume corrections 
to the asymptotic value differ for both ensembles. 
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Finally, using the results of the probability functions we have derived the asymptotic behavior of 
the charged particle moments and established the differences in the grand canonical and canonical 
formulation. We have also applied these results to find the thermodynamic limit of a class of semi- 
intensive quantities. It was shown that in systems with exact and average charge conservation 
such quantities should naturally converge to different values in the thermodynamic limit. This is 
because the behavior of the semi-intensive quantities in the near vicinity to the thermodynamic 
limit are determined by the subleading, finite volume, corrections to the probability distributions 
which are specific to a given statistical ensemble. 

Is important that first moments are the same in the canonical and grand canonical ensemble. 
This means that particle yields in heavy ion collision and equation of state of dense hadronic 
medium arc insensitive to the statistical ensemble in the thermodynamic limit. This is not the 
case, however, for fluctuations and higher moments. Finite volume effect are more and more 
relevant for higher moments. Such a situation appears when comparing the statistical model with 
lattice gauge theory results obtained on a small lattice. 



[1] For a review see, e.g. P. Braun-Munzinger, K. Redlich, J. Stachel: Quark Gluon Plasma 3 eds. 

R.C. Hwa, X.N. Wang (World Scientific, Singapore 2004) 491-599; A. Andronic, P. Braun-Munzinger: 

Lect. Notes Phys. 652 35 (2004) 
[2] J. Cleymans, K. Redlich, L. Turko: Phys. Rev. C 71 047902 (2005) 
[3] J. Cleymans, K. Redlich, L. Turko: J. Phys. G 31 1421 (2005) 

[4] V.V. Begun, M. Gazdzicki, M.I. Gorenstein, O.S. Zozulya: Phys. Rev. C 70 034901 (2004); V.V. Be- 
gun, M.I. Gorenstein, A.P. Kostyuk, O.S. Zozulya: Phys. Rev. C 71 054904 (2005); V.V. Begun, 
M.I. Gorenstein, O.S. Zozulya: Phys. Rev. C 72 014902 (2005); A. Keranen, F. Becattini, V.V. Begun, 
M.I. Gorenstein, O.S. Zozulya: J. Phys. G 31 S1095 (2005) 
[5] K. Redlich and L. Turko: Z. Phys. C5 201 (1980); L. Turko: Phys. Lett. B 104 153 (1981) 
[6] CM. Ko, V. Koch, Z. Lin, K. Redlich, M. Stephanov, X.N. Wang: Phys. Rev. Lett. 86 5438 (2001); 
Z. Lin and CM. Ko: Phys. Rev. C 64 041901 (2001); S. Jeon, V. Koch, K. Redlich, X.N. Wang: 
Nucl. Phys. A 697 546 (2002); K. Redlich, V. Koch and A. Tounsi: Nucl. Phys. A 702 326 (2002); 
K. Redlich, J. Cleymans, H. Oeschler, A. Tounsi: Acta Phys. Polon. B 33 1609 (2002). 



